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Abstract 

We study dynamics of the shearless stratified turbulent flows. Using the method 
of differential constraints we find a class of explicit solutions to the problem under 
consideration and establish that the differential constraint obtained coincides with 
the well-known Zeman-Lumley model for stratified flows. 

1 Introduction 

The results of this article continue our earlier studies of the problem of interaction and mi- 
xing between two semi-infinite turbulent flow fields of different scales given in [Q, ||] wherein 
it was proposed a concept based on the method of differential constraints [23, |H| for exa- 



mining the closure procedure for momentum equations in Parametric Turbulent Models. 
The key idea of this approach can be formulated shortly as follows: the algebraic expres- 
sions for the n-order moments of statistical characteristics of turbulent flow are determined 
as the equations of invariant sets (manifolds) of the corresponding differential equations. 
As an illustrative example we have shown how this concept can be applied for finding 
a selfsimilar solution to the above-mentioned problem in the case of nonstratified flow. 

The aim of the present article is to extend the proposed method to the general case 
including stable and unstable stratification of the flow. At first, we find a differential 
constraint generated by the model under consideration and then we construct a reduction 
that enables us to rewrite that model in a more simple form on the invariant set obtained 
both for stable and unstable stratifications. This makes it possible to find explicit solutions. 
The analysis of those solutions shows that the influence of stratification on statistical 
characteristics of turbulent flow is essential. The time scales of turbulence for the all 
cases are founded in explicit forms and we show that their behaviors are different. As an 
application to the Theory of Parametric Turbulent Models, we indicate that the differential 
constraint obtained coincides with the well-known Zeman-Lumley model ||24||. 
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The motivation for the study is that, in many practical situations, a turbulent flow 
is embedded in a surrounding field of turbulence of different intensity. The shearless 
turbulent flows appears in decaying grid turbulence in which the mean velocity is constant 
throughout. Those flows include both the homogeneously-damped grid turbulence and 
shearless turbulent mixing layer which is formed beyond a composite grid. 



2 The governing equations 

The directions of turbulent-flow investigations based on constructing closed equations 
of turbulent transport are used intensely in modelling. The approach based on two- 
parametric models of turbulence gained wide application, as well as one based on second- 
order closure models (for example, the K — e model), which are effective from the compu- 
tational viewpoint and yields results whose accuracy is sufficient for practical applications. 
That models are based on parametrization of third-order moments of the gradient type. 
However, the use of that models for description of turbulent transport in stratified flows 



gives a qualitatively incorrect result (in some cases) see, for example [16]. The anisotropic 



character of the buoyancy effect on the structure of turbulence is manifested in appear- 



ance of the long wavelength spectrum of turbulent oscillations [15]. This spectral range 
corresponds to large-scale vortex structures containing the main portion of turbulent en- 
ergy. According to experimental and theoretical studies, the following large-scale vortex 
structure is formed in stratified flows: turbulent spots in the case of stable stratification 
and coherent structures in the case of unstable stratification which are mainly responsible 
for turbulent transport. The effect of intermediacy and asymmetry of vertical turbulent 
transport caused by the influence of large-scale vortex structures makes the probability 
distributions of turbulent fluctuations significantly non-Gaussian. The turbulent struc- 
ture in these flows is usually described by third-order closure models where the triple 
correlations (asymmetry) are calculated from differential transport equations. There is 
considerable number of references connected with the use of third-order turbulent mo- 
dels. However, as a rule these models employ Millionshchikov's quasinormality hypothesis 
for the parametrization of diffusion processes in equations for triple correlations and, ac- 
cording to this hypothesis, all cumulants of fourth- and higher-order may be negligibly 
small in comparison with the corresponding correlation functions. As a consequence, in 
the former case the triple-correlation equations are of the first-order without a dump- 
ing mechanism for triple correlations that leads to physically contradictory results flq| . 
The approach proposed in [1C, [Tl|] allows us to overcome this obstacle; the technique also 
includes a physically reasonable way for constructing approximate algebraic parametriza- 
tions of higher moments. Using these observations, we present a third-order model of 
turbulence to describe correctly the shearless turbulent mixing layer in the framework 
stated above. 

To determine the values of the horizontal components eh of the turbulent kinetic energy 
e = e^, + l/2(w 2 ) and the rate of dissipation e and the one-point velocity correlation (w 2 ) 
of the second-order, we make use of the differential equations: 

de h _ d(e h w) c x \ _2 1_2 

at ~ dz t eh 3 3 e ' 
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de d 
dt dz 
d(w 2 ) d(w 



c d r(w } 



2 > de 



dz 



T 



dt 



dz 



+ 2Pg(w9) 



£l 

T 



-E 



-€. 



Here E, r = E/e are the kinetic energy and the time scale of turbulence respectively. 
The volumetric expansion coefficient is (3 = 1/©; and 9 are the mean and variance 
potential temperatures respectively. The constants involved in the model with the lower 
case letters are denoted by c**. The system of equations is nonclosed and we complete it 
by the transport equation for the triple correlation of the vertical velocity fluctuation: 



d{w 3 
~~dt 



, 9. d(w 2 ) „ , o 



dC 

dz 



+ 3f3g(w z 9) - c 2 



(w 3 ) 



where C = (w 3 ) — 3(w 2 ) 2 is the cumulant of velocity fluctuations. To obtain a closed 
model of turbulent transport that does not imply equality to zero of the fourth-order cu- 
mulants, the closure procedure is performed at the level of the fifth moments, i.e. it is 
assumed that the fifth-order cumulants are equal to zero. Thus, to obtain the distribu- 
tion of the third-order cumulants (correlations), the latter are calculated from differential 
transport equations: the fourth-order cumulants are determined approximately (from al- 
gebraic expressions), and the fifth-order cumulants are assumed to be equal to zero, since 



their contribution is negligibly small. Results \12\ of numerical simulation of vertical tur- 
bulent transport in a convective boundary layer confirm the validity of this approach. To 
express the fourth-order cumulant C, we can use the equality 



c = - T - 

C3 



, , v d(w 2 ) , 9v d(w 3 ) 



dz 



dz 



The triple correlation (w 2 9), the vertical heat flux (w9) and the temperature dispersion (9) 



are approximated by the following algebraic relationships [11]: 



( w 2 9) = - — 
c 4 



dp) 

( w 3 )—-2(3g(w9 2 



t , o.de tn 2 



dz 

dz 
dQ 



(w6 



if 



C$1 

— (w6) 
Cg 1 r oz 



/3gce 1 
tN 2 
(3gr 



<- 2 >, 
(w9) 



N 



C5 

de 

dz 



d(9 2 



dz 



where r = t/tq, tq is the time scale of potential temperature variance, N is the Brunt- 
Vaisala frequency. On using the balance approximation between exchange mechanism and 
dissipation, the equation for the horizontal component of the turbulent kinetic energy 
may be simplified and may be written as follows: 



-ci 



eh 



eh + 



(w 2 ) 



Hence, 



eh 



ci 



1 



ci + 2 



E 



3ci 
2(ci + 2) 



eh + 



(w 2 ) 



3ci (w 2 ) 
2(ci + 2) e 
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As a result of the above simplifications, the model of the shear less mixing turbulence layer 
of third-order is represented by the following system of differential equations 



d(w 3 
~~dt 



d_ 

Oz 



T 
C3 



d(w 2 ) _ d(w 3 ) ci (w 2 ) 
dt ~ dz ci + 2 t* ' 
(dw 3 



6(w 3 )^^- + A(w 2 
dz 

de d 
dt dz 



Oz 



3K) 



2 .d(w 2 ^ 



Oz 



C2- 



(w 3 ) 



C d T(w 2 ) 



de 
Oz 



Here 



T 



6(ci+2) 
( Cl +4)c ei r 2 Af 2 



1 + 



1 + 



C2C4T 



c E1 (ci+C2) 
Ce2 ( Cl + l)c ei r2V2 



1 



It follows from the formula for C |IT| that the contribution of the second term in the 
algebraic model for the cumulant C is essential. Also, we can assume that (see |ffif| ) 



r* = t** = t*** 



Thus the governing equations are: 

d(w 2 ) d{w 3 

d(w 3 ) 



dt 



dz 



1 + ^T 2 N 2 



a(l + at 2 N T 



(w 2 ) 



dt 



d_ 

dz 



kt(w ) 

de d 
dt dz 



d(w 3 
dz 



St(w 2 



-3(w 2 ) 

de 
dz 



d(w 2 



7 (1 + at z N 



1 v <™ 3 > 



dz 

g(l + at 2 N 2 ) -, 



(2.1) 
(2.2) 

(2.3) 



where f = (w 2 )/e and a = 2/3, k = 601/03(01+ 2), 7 = 2c2(ci + l)/3ci, 5 = 3ciq/2(ci + 2) 
q = 2c e2 (ci + 2)/3ci, a = cfvr/18(ci + 2) 2 . 
In addition, we indicate the equation for f 



dt 
dt 



{w 2 ) 



d(w 3 } 2\d 2 ( w ' 2 ) f 9{w 



dz 



+ 5t{w* 



dz" 2 



+ St 



2\\ 2 



dz 



2 - d2f , ^ x ~ dfd ( w 



<i(^) (io 2 ) + (f>-a)(l+< I f 2 A' 2 ) (2.4) 



<9~ 



which can be obtained from fl2.1|), ( |2.3D . This equation will be crucial for the further study 
of properties of system d2.lD- (|2.3| ). 



3 Invariant sets 

In this section we show that system ( |2.1| )-( plj|) admits an invariant set (manifold) of the 
form 

D = {(w 2 ), (w 3 ),t : H\{w 2 ), (w 3 ),t) = (w 3 ) + 5t(w 2 )(w 2 ) z = 0} . (3.1) 
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A possible application of the existence of an invariant set (manifold) generated by the 
differential equation for the triple correlation of the vertical velocity fluctuations consists 
in constructing explicit solutions to the model obtained. Here "explicit" means solutions 
which can be found by means of ordinary differential equations (or algebraic systems). 

Let us briefly present the special terminology of Symmetry Analysis (see |9], [ll| for 
more details). 

Consider a system of evolution equations J- 



u\ — J- 1 (t, Xi, . . . , x n , u 1 , . . . , u\, . . . \ 



where i = 1, . . . , m, u\ = d x u k /dx± 1 ■ ■ ■ dx Xn . 
A set (manifold) 7i given by equations 

hi (t, x\, . . . , x n , . . . , u , ... ,u , . . . , 

is said to be the invariant set (manifold) of system T if 

vArf)^ =o, 



dt 4—1 du 1 4-f du\ ' 

t=i i=i 

where a = (a 1 ,..., a n ), D a = D%\ ■■■D%™. 

The invariant condition can be written in the following equivalent form 

D t {hi) \ [^ \ [n]o = 0. 

Here [J^q is oo-prolongation (see |l^]) of T with respect to xi,...,x n . The set [H]o is 
determined by analogy. 

To show that D is invariant under the flow generated by equation (|2.2|) , we prove 
that the operator TL l {(w 2 ), (ui 3 ),f) = (w 3 ) + 8t(w 2 )(w 2 ) z preservs the sign on the set of 
sufficiently smooth solutions to system (|2 . 1| ) — (|2 . 3|) . 

We introduce into consideration the following sets 

D + = {p, q,we C°°(R) : H\p, q, w) > 0} , (3.2) 
D~ = {p, q,we C°°{R) : H\p, q, w) < 0} . (3.3) 

Definition 1. Operator TC 1 is said to be sign- invariant of system ( ^.l|) - (|2.3| ) if 

((w 2 (;t )),(w 3 (;t )},f(;to)) G D+(D~) 
=> {(W 2 (;t)),(w 3 (;t)),f(;t)) £ D+(D-), t > t . 

Remark 1. For quasilinear parabolic equations derivation of suitable sign- invariants plays 
an important role for existence, uniqueness, regularity, existence of explicit solutions and 
other problems (see, for example, [|l4|, §). 
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It is convenient to determine the invariant set D under the flow generated by equa- 
tion ( [2.2D via the sign-invariant sets D + and D~ as their intersection. The invariance 
condition takes the form 



d_ 

at 



n 



o. 



As the first result that uses the above notion we have: 

Theorem 1. Let {((w 2 ), (w 3 ), e)} be a set of sufficiently smooth solutions of (2.1)-$2j\) 
such that the functions f = (w 2 )/e satisfy the relations 



df 

7T = ' 

dz 



| = (2a- 7 )(l + af 2 ^) + ^. 



(3.4) 



Assume that k = 5. Then operator Tt 1 is a sign-invariant of system ([IT 
Proof. Calculating the time derivative, we obtain 



d(w 3 ) r d L . n. d(w 3 ) 



dt 



dz 



(3.5) 



Using equation (|2.2j) and the assumption that df/dz = 0, we can rewrite (3.5) as follows 



d_ 

dt 



n 



d_ 

dz 



d(w 3 
dz 



3(w 2 ) 



d(w 2 



dz 

d(w 2 ) d(w 2 
dt dz 



7 M (l + ar 2 N 2 ) 



+ Sf{w 



dtdz 



Replacing the derivatives d(w 2 )/dt and d 2 (w 2 } /dtdz by their representations from equa- 
tion (p.l|) (taking into account that df/dz = 0), we have 



d_ 

dt 



n 1 



d(w 2 ) d(w 3 ) 

Kf \ \ + KT{W 2 ) 



d 2 {w 3 



dz 



dz 



df d(w 2 ) . 2 . 
+ 5 dt^dz 2{w) 

+ Sf(w 2 ) 



Sf 



dz 2 
d(w 3 



3{w 2 ) 



d{w 2 
dz 



7 ^ (i + a f^N 2 ) 



d 2 (w 



dz T 
3v ad(w 2 ) 



(l + af 2 N 2 ) 



d{w 2 
dz 



dz 2 



f dz 



af 2 N 2 ) 



This formula can be written in the form 



d_ 

dt 



KT- 



d(w 2 ) d(w 3 ) 
dz dz 



+ Kf(w 2 ) 



d 2 (w 3 
dz 2 



3<tO 



2 .d(w 2 ^ 



In 1 (1 + at 2 N 2 ) + 7 <5(1 + at 2 N 2 )(w 2 )^^ + S 



d(w 2 



St 



d{w 2 ) d{w 3 



dz dz 



2Sa (1 + ar 2 N 2 ) 



d{w 2 
dz 



dz 

-(w 2 )-5f 



dz 

.df d(w 2 
dt dz 
d(w 2 } d 2 (w 3 
dz dz 2 



(w 2 ) 



As a result, we obtain 



at dz dz 



/ CN-/ 2^ 2 (W 3 ) . r,,d(W 2 ) 

+ (k-5)t(w 2 )^-^ + (w 2 )-^ 



2a5 (1 + af 2 N 2 ) + jS (l + af 2 A^ 2 ) - 3) - 1 (l + at 2 N 2 ) H 1 



dz 2 x 1 dz 

T 
t 



5 Tt 



(3.6) 
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It follows from the equality k = 8 and (|3.4D that 



a 



1LU 1 = - (l + af 2 N 2 ) H 1 . 

ut T ' 



Therefore 



n 1 



t=ti 



n 1 



t=to 



exp 



to 



+ of N 2 ds 



This completes the proof of the theorem. 



Theorem 2. Let a(2a — 7) = a(g — a), | + 2a 



7 = £ — a and k = 8. Then system 



(2.1)-(2.S) admits the invariant set D and its reduction on the set D is of the form: 



(w 2 ) = fe, 



( w 3) = -8f(w 2 ) 



d(w 2 



de 

at 



d_ 

Dz 



ot(w ) — 
dz 



dz 

g(l + af 2 N 2 ) % 



where the function f(z,t) = f(t) solves the ordinary differential equation 



df . 
~dt ={e 



a) 



'I + af 2 N 2 ) 



(a version of equation (EOI) on the set D) 



(3.7) 
(3.8) 

(3.9) 
(3.10) 



Proof. Observe that, according to Theorem |l| the invariant set exists if the equation for f 
admits a class of solutions which satisfy conditions ( |3.4| ). To find this class of solutions, we 
consider equation (|2.4[) for f . As we have indicated above, this equation is a consequence 
of equations ([O]), (pT3|). More precisely, calculating the time derivative for f , we obtain 



Of 
dt 



ld{w 2 ) {w 2 )de 



dt 



8{w 2 )f 



d 2 i 



dt 



1 d(w 3 
e dz 



a(l + af 2 N 2 ) 



8f 



,d 2 (w 2 



dz 2 



dz 2 {w 2 ) \ dz 
Thus the equation for f is of the form 



d{w ^ +0 (l + at 2 N 2 ) + 25 dfd{w2) 



dz dz 



df 
dt ~ 

+ 5(w 2 )f 



(w 2 ) 



d{w 3 ) , ^ t ^ 2 \d 2 {w 2 ) 



d 2 ^ 



+ 2<5f 



dz 

df d(w 2 
dz dz 



+ 8f{w z 



dz 2 



df 



+ 8f 



d(w 2 
dz 



dz 



{w 2 ) + (g-a) (l + af 2 N 2 ) . 



Obviously, it is sufficient to check the conditions of Theorem |T| only for ((w 2 ), (w 3 ),f) G D. 
It is clear that equation (|2.4|) on the set D can be rewritten in the form 



df 
~dt 



(g-a) (l + af 2 N 2 ) 



for f(z,t) = f(t) (the expression in square brackets equals zero). Using the equalities 
a(2a — 7) = a(g — a), | + 2a — 7 = g — a, we conclude that the above-mentioned equation 
and equation defined by ( |3.4[ ) coincide with each other. The proof is completed by simple 
checking that the function (w 2 (z,t)} = f(t)e(z,t) satisfies identically (2.1) on the set D, 
where f solves the equations ( |3.10| ). ■ 
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Theorem [l] is of a special interest in view of its application to Turbulent Models. 



Corollary 1. The equation TL 1 {{w 2 ), (u> 3 ),f) = that defines an invariant set of (2.1) 



(2.5) coincides with the algebraic triple correlation model or the Zeman-Lumley model t 2_4] 



In other words, the algebraic expression represents the equation of an invariant set 
(manifold) generated by the differential equation for the triple correlation. 

4 Solutions on invariant sets 

Theorem||| enables us to reduce (|2.1[ )- (|2.3| ) to the algebraic differential expressions (jOI)- 
( 3.1 0| ) which can be easier analyzed. Using the obtained reduction, we construct explicit 



solutions to system fl2.1|) -([2.3|) for iV 2 ^ 0. In the case of nonstratified flow, i.e. for 
iV 2 = 0, it was proven in |5[ that system (2.1)-( [2^ ) admits a parametric group of scale 



transformation that enables us to find a selfsimilar solution of the form 

HO , 2n m > 3 N 



a (t + t ) 3 ^' OjI ( t + tQ )2 M ' (t + t ) 3 ^ 

£ = ^-^, L = \(t+t ) U , Z C = \ L + \ 1 , t >0, 

where A, Aj are model constants, to is a parameter and u = 1 — fx. In |3], |5[ we studied the 
existence of selfsimilar solutions and presented in detail the qualitative properties of the 
solution obtained. 

A direct calculation shows that for stratified flows there are no selfsimilar solutions 
similar to e a , (ii> 2 ), {w 3 ). Nevertheless we find a class of explicit solutions (^1|)-( |2^) by 
using the fact that system (fP|)-(|!0|) is equivalent to algebraic differential expressions 
-( |3.10| ) under the above-mentioned hypotheses on the parameters of the model. 
Let us rewrite ( |3.10| ) in the form 

jL = BN 2 f 2 (t) + D, (4.1) 



where B = a(2a— 7), D = 3/5+2a—j. We note that the form of solutions to ( |4.lD depends 
on the signs of quantities B, D and iV 2 . Positiveness (negativeness) of iV 2 corresponds 
to the case of stable (unstable) stratification. The signs B and D are determined by 
values of coefficients c** which are experimentally known numbers. The calculated values 
of numbers B and D show that B < 0, D > 0. 

4.1 Stable stratification 

Integrating (^[l]) and denoting by 



D 



-BN 2 

we can obtain the positive solution f s = f s (t) which is defined by the formula 



f s (t) = Atanh Uy/{D/ - B)BNt + C , 



152 



V N Grebenev and B B Ilyushin 



where Co is a constant. Given the initial data f s (to), the constant Co can be easily 
determined. We note that function f s (t) has a horizontal asymptote; more exactly, r s (t) — ► 
^4 as t — * oo. 

Our main aim is to find a solution to fl3.9|) . The initial condition for equation fl3.9|) is 
determined by the physical model. We have 



e (z) = e(z,Q) 



e_, if z < 0, 
e_|_, if 2 > 0, 



where e_, e+ are the positive numbers such that e_ ^ e + and we assume that e_ < e+. 
Set 



9(t) 



and 



e(z, 



u s {z, V) exp 



tj)(p)dp J , where e(z, 0) = e(z, i). 



The function 0(t) maps [0, +oo) onto [0, +oo) and for u s we have: 



du s 
39 



D_ 

dz 



5u s 



du s 
dz 



where 9 = J exp ^— J if)(p)dp ) 
u s (z,0) = e+ if z > 0. 



u s (z, 0) = e_ if z < 0, 
It is easy to check that 6 : [0, +oo) — > [0, 9q), where 9q < A/aN 2 and 



(4.2) 
(4.3) 



exp 



i[i(p)dp 







as 6* — > +oo. 



In studying the Cauchy problem for equation (3.9), we base our analysis on investigation 
of transformed problem (4.2), (4.3) for a finite time interval [0, #o)- Equation ( f4.2j ) is 
usually called the porous medium equation. It is well-known that solutions to (|4.2|), (f4.3| ) 
are unique and invariant under a parametric group of the scale transformation (see pOfl ). 
Therefore u s is a selfsimilar solution which can be represented in the form 



u s (z,e) = u s (i), | 



Z — Z r . 



z c 0) = A O V20, 



where Ao is a model constant, and then fl4.2| ), (4.3) is rewritten as 

2 



d£ 2 V d£ 



. * , s du s 

(e + A )-/ = o, 

eft; 



(4.4) 
(4.5) 



u s (-oo)=e-, u s (+oo) = e + . 
Here z c (#) is the so-called central line. Setting u s (£) = d^/dC, where £ is a new variable 



we obtain the following boundary value problem for the Blasius equation [15] 
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Equation (|4.6| ) arises in the context of the evolution of boundary layers in an incompressible 
fluid along a surface |j"H| . As a result, we have that there exists an one-parametric family 
of solutions to (4.6), (4.7). This implies existence of a solution to ( [4.4[ ), (|4.5|), In Q we 



proven this result directly for a problem of type (4.4), (4.5) and noted that this solution 



is essentially different from the well-known Barenblatt's solution 

Remark 2. For the first time, the class of positive solutions to the porous medium 



equation was introduce in p7|. 



Thus we arrive at the following 

Lemma 1. For any positive finite e_ and e+ (e_ < e+ ) there exists a unique positive 
solution u s to &4-?V (respectively (J^i), such that junction u s is increasing 



over (— do, +00J; moreover u s has convex and concave profiles for £ < and £ > 

Remark 3. The test configurations of profiles of the spectral flux e obtained by numerical 
and experimental methods (see, for example jl], |2l|) coincides qualitatively with profiles 
of u s . 

Once we have determined u s , we can find (w 2 ). The function (w 2 ) is defined from the 
relation 

(w 2 (z,6)) =U{z,0) =f s u s (z,6)exp (- [ i>(p)dp) , {w 2 (z,t)) = (w 2 (z,6)). 



\ Jo / 

For the triple correlation (w 3 ) we obtain 

(w 3 (z,e)) = -5f s {w 2 (z,e)) d{w2 ^ e)) , {w 3 (z,t)) = {w\z,8)). (4i 
4.2 Unstable stratification 



Let us now find a solution to ( ^4.1[ ) for A^ 2 < 0. In this case we arrive at the following 
solution to equation (|4.1|) 



f ns = ,4 tan \J (DBN 2 )t + Ci 



where G\ is determined by the initial data f (0). As before, we introduce into consideration 
the new time variable 

e = e(t)= [ T 2 s (p)d P 

Jo 

and note that 6 — > +00 as VDBNH + d vr/2. Then we obtain 

u ns (z,0) = e(z,6)exp (J ip(p)dp^J , 

m = 1 + a f(9) N ^ <o)=f 2 ns (e-Ht)), 
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and 



du r , 



06 



dz 



c du ns 
ou ns — — 
oz 



where 



cxp 



ip(p)dp I d£. 



Here 9 -> +oo as y/(DBN 2 )t + d -» vr/2 and exp (- £ tp(s)ds) is an increasing bounded 

function on the interval [0, +oo) that coincides with experimental observations about in- 
creasing the spectral flux e in the case of unstable stratification of the flow. Further 
analysis goes along the same lines just as for N 2 > 0. Combining Lemma [l] with the above 
arguments, we claim the following 



Theorem 3. Let k = 5 and 

a(2a — 7) = a(g — a), 



+ 2a — 7 = g — a. 



Then there exists a solution to system (2.1)-(2.S) of the following form 



e[z, 



u s (z, 0) exp 



(w 2 (z,t)) = f s e(z,t), 



and 



i(z, 0) = u ns (z, 9) exp 
(w 2 (z,t)) = f ns e(z,t), 



ip(s)ds 
(w 3 (z,t)) 

ijj(s)ds 
(w*(z,t)) 



0(1) = j t 2 (p)dp, e(z,t)=i(z,9), 
d{w 2 (z,t)) 



-5f s (w 2 (z,t)) 



dz 



for N z > 0, 



9 = 0(t) = 

-dT ns (w 2 (z,t)) 



Tns(p)dp, e{z,t) = e{z,9), 
d{w 2 (z,t)) 



dz 



for N z < 0, 



where u s (u ns ) is a selfsimilar solution to ( \4-4 )> l U-'V - 



5 Conclusions 

As a rule, the parametric models of turbulence represent the transport equations for the 
second-order moments which are completed by closure relations for the higher-order mo- 
ments and dissipation tensor, and in many cases these closure relations are given in the 
so-called isotropic form. The assumption about the relaxation character of turbulence un- 
der its evolution to the equilibrium state (the homogeneous isotropic state with Gaussian 
distribution of turbulent fluctuations) is the base for such simplification. However, these 
models describe adequately certain statistical structures of investigated flows even if tur- 



bulence is characterized by anisotropic effects. In general, conditions (3.4) (which assume 
"uniformity" of the turbulent scales) can be used in the case of an equilibrium state of 
turbulent flows. Nevertheless, the results obtained provide a correct modelling for turbu- 



lent flows: the algebraic model for triple correlations (|3.8|) was used in [24] for studying 



turbulent structures in the convective boundary layer and moreover, it was testified for 
distinct turbulent flows (see [fTcR ). 



We conclude with some observations and comments. A version of formula (3.8) coincides 



with the well-known Hanjalic-Launder model M in the case of nonstratified flow Mj. It is 



Invariant Sets and Explicit Solutions 



155 



possible to get results similar to the results of this article for shear flows in the problem 
of plane turbulent wake. Maybe the most interesting and novel part of the criterion of 
invariance for the so-called locally equilibrium approximations (see [Q) in the problem 
of plane turbulent wake consists in showing that the Poisson bracket {U, e} equals zero. 
Here U is the velocity excess and e is the turbulent kinetic energy. An example of such 
an equality was obtained in for a selfsimilar solution. 



Acknowledgements 

This research was partially supported by Integration Project SD RAS (grant no. 2000- 
01), RFBR (grant no. 01-01-00783). This work was supported by INTAS (proposal no. 
97-2022). 



References 

Alexeenko N V, Bukreev V I and Kostomakha V A, The Shearless Interaction of Two Isotropic 
Turbulent Fields, J. Appl. Median. Tech. Phys. 1(149) (1985), 57-62. 

Barenblatt G I, Scaling, Selfsimilar and Intermediate Asymptotics, Cambridge Texts in Ap- 
plied Mathematics, Vol.14, 1996. 

Galaktionov V A, Quasilinear Heat Equations with First-Order Sign-Invariants and New 
Explicit Solutions, Nonl. Analysis TMA 23(12) (1994), 1595-1621. 

Grebenev V N, Ilyushin B B and Shokin Yu I, The Use of Differential Constraints for Ana- 
lyzing Turbulence Models, J. Nonl. Sci. Numer. Simulation 1(4) (2000), 305-317. 

Grebenev V N and Ilyushin B B, Application of Differential Constraints for Analyzing Tur- 
bulence Models, Dokl. Russian Acad. Sci. Mech. 374(6) (2000), 371-375 (in Russian). 

Grebenev V N, On Dynamics of Far Plane Turbulent Wake, Sib. J. Ind. Math. 1(1) (1998), 
97-103 (in Russian). 

Hanjalic K and Launder B E, A Reynolds Stress Model of Turbulence and Its Application to 
Thin Shear Flows, J. Fluid Mech. 52 (1972), 609-638. 

Hanjalic K and Launder B E, Fully Developed Asymetric Flow in a Plane Channel, J. Fluid 
Mech. 51 (1972), 301-335. 

Ibragimov N H, Transformation Groups Applied to Mathematical Physics, Reidel, 1985. 

Ilyushin B B, Higher-Moment Diffusion in Stable Stratification, in Closure Strategies for Tur- 
bulent and Transition Flows, Editors: Launder B E and Sandham N D, Cambridge University 
Press, in press. 

Ilyushin B B, Model of Fourth-Order Cumulants for Prediction of Turbulent Transport by 
Large-Scale Vortex Structure, J. Appl. Mechan. Tech. Phys. 40(5) (1999), 871-876. 

Ilyushin B B and Kurbatskii A F, Numerical Simulation of the Shearless Turbulence Mix- 
ing Layer, Izv. Acad. Nauk USSR (Siberian Branch), Ser. Tech. Nauk 3 (1990), 62-68 (in 
Russian). 

Kaptsov O V, B-Determining Equations: Applications to Nonlinear Partial Differential Equa- 
tions, Euro. J. Appl. Math. 6 (1995), 265-286. 



156 



V N Grebenev and B B Ilyushin 



[14] Kalashnikov A S, Some Problems of the Qualitative Theory of Non-Linear Degenerate Second- 
Order Parabolic Equations, Usp. Mat. N auk 42 (1987), 135-176 (in Russian); English trans- 
lation: Russ. Math. Survs., 42 (1987), 169-222. 

[15] Kochin N E, Kibel I A and Rozc N V, Theoretical Hydromechanics, Part II, Moscow, FML, 
1963 (in Russian). 

[16] Lamb R G, Diffusion in Convective Boundary Layer, in Atmospheric Turbulence and Air 
Pollution Modelling, F.T.M., Editors: Nieuwstadt and H. van Dop, Reidel, Boston, MA, 
1982. 

[17] Polubarinova-Kochina P Ya, On a Nonlinear Partial Differential Equation Appearing in the 
Filtration Theory, Dokl. Acad. Nauk USSR LXIII(6) (1948), 623-626 (in Russian). 

[18] Monin A S and Yaglom A M, Statistical Hydromechanics, Moscow, Nauka, 1967 (in Russian). 

[19] Sidorov A F, Shapeev V P and Yanenko N N, The Differential Constraints Method and Its 
Application to Gas Dynamics, Nauka, Novosibirsk, 1984 (in Russian). 

[20] Vazques J, The Interfaces of One-Dimensional Flows in Porous Media, Trans. AMS 285(2) 
(1984), 717-737. 

[21] Veeravalli S and Warhaft Z, The Shearless Turbulence Mixing Layer, J. Fluid Mech. 207 
(1989), 191-229. 

[22] Wcinstok J, A Theory of Turbulence Transport, J. Fluid Mech. 202 (1989), 319-338. 
[23] Yanenko N N, Selected Works, Nauka, Moscow, 1991 (in Russian). 

[24] Zeman O and Lumlcy J L, Modeling Buoyancy Driven Mixed Layer, J. Athmos. Asci. 33 
(1976), 1974-1988. 



